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Ultracold neutrons (UCNs) were produced in superfluid helium using the PF1B cold neutron beam 
facility at the Institut Laue-Langevin. A 4-liter beryllium-coated converter volume with a mechanical 
valve and window-less stainless steel extraction system were used to accumulate and guide UCNs to a 
detector at room temperature. At a converter temperature of 1.08 K the total storage time constant 
in the vessel was (20.3± 1.2) s and the number of UCNs counted after accumulated was 91,700±300. 

From this, we derive a volumetric UCN production rate of (6.9 ± 1.7)cm _3 s _1 , which includes a 
correction for losses in the converter during UCN extraction caused by the short storage time, but 
not accounting for UCN transport and detection efficiencies. The up-scattering rate of UCNs due 
to excitations in the superfluid was studied by scanning the temperature between 1.2-2.4 K. Using 
the temperature-dependent UCN production rate calculated from inelastic neutron scattering data, 
the only UCN up-scattering process found to occur was from two-phonon scattering. Our analysis 
for T < 1.95 K rules out the contributions from roton-phonon scattering to < 29% (95% C.I.) and 
from 1-phonon absorption to < 47% (95% C.I.) of their predicted levels. 


I. INTRODUCTION 

Ultracold neutrons (UCNs) are neutrons with kinetic 
energies less than the neutron optical potential U op t of 
well-chosen materials; for instance, beryllium has U^ e = 
252 neV. They can reflect from material surfaces at all 
incident angles, allowing them to be stored in a vessel 
and studied for times up to a few times the neutron 
lifetime (880.3 ± 1.1) s pQ. Due to their low velocities 
(v ~ 5ms _1 ) and the long times they can spend in media 
or fields chosen by an experimenter, UCNs have become a 
valuable tool for high-precision studies of the fundamen¬ 
tal properties of the neutron and its interactions. These 
studies have wide-ranging applications in nuclear physics, 
particle physics, and cosmology |2H5]. For instance, they 
are currently used in searches for the permanent electric 
dipole moment of the neutron EHH, measurements of 
the neutron lifetime [T5H25] and /3-decay correlation pa¬ 
rameters [241 : '25 \, and quests for dark matter candidates 
PHHliHj, axion-like particles 1201-1.1?). and Lorentz invari¬ 
ance violations [551 . 

These experiments would benefit from an increased 
density of UCNs, which has motivated the development 
of “next-generation” UCN sources 1541155) . They con¬ 
vert cold neutrons (CNs) to the UCN energy range by 
allowing them to generate excitations in solid deuterium 
or superfluid 4 He [44] [45]. The latter has the advan¬ 
tage that 4 He has a zero neutron absorption cross-section 
and, if the converter is kept at sufficiently low tempera¬ 
tures (typically < 1 K), thermal up-scattering of UCNs is 
sufficiently suppressed. This allows the produced UCNs 
to survive in the converter for times dominated by wall 
losses of the vessel, typically ~ 100 s MM- Since UCNs 
are not in thermal equilibrium with the converter, densi¬ 
ties can be greater than those possible from direct ther¬ 
mal moderation, hence the name “super-thermal” pro¬ 


duction. 

The need for more intense UCN sources has motivated 
the present work of developing a super-thermal helium 
source with a mechanical UCN valve separating the con¬ 
verter vessel from the extraction guides. This allows 
UCNs to be accumulated and then used to fill the volume 
of an external experiment when it is required. A unique 
feature of our source is the ability to extract UCNs from 
the converter using a vertical extraction system so that 
windows, which can cause a loss of UCNs, are avoided 

m- 

Locating the UCN converter volume and its associated 
cryogenics at the end of a CN beam guide, far from the 
initial source of neutrons, typically from spallation or 
fission processes, reduces the needed cooling power and 
activation of the source. Furthermore, the source can be 
adapted for each experiment in order to optimize UCN 
delivery; for instance, the length, volume and coating of 
the UCN transport guides and converter vessel can be 
altered. 

A successful realization of such a superfluid helium 
source, which is dedicated to an experiment investigat¬ 
ing the neutron quantum levels in Earth’s gravitational 
field [49], is described in Refs. [41] 02]. The still large 
room for improvement has motivated the development 
of a second UCN source at the Institut Laue-Langevin 
(ILL), Grenoble, France. The first experiments testing 
and characterizing this source using the PF1B CN beam 
facility are presented in this paper. 

Inelastic scattering of neutrons [S3 El] has been used 
to reveal the fundamental excitations in superfluid he¬ 
lium [52] [55]. These studies can also be performed by 
measuring the up-scattering of UCNs caused by interac¬ 
tions between a neutron at rest and excitations in super¬ 
fluid helium at different temperatures T. Using UCNs, 
one is more sensitive to weak scattering processes due to 
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the intrinsic 1/v dependence of the cross section, as well 
as the potential hundreds of meters path length a UCN 
can spend in the medium under study [53j . 

The contributions to the up-scattering rate constant 
r^p 1 (T) for T < 1 K have been determined [35] using the 
Landau-Khalatnikov Hamiltonian [55] [56] and the three- 
phonon interaction described in Ref. m to be 

Rip (-^) T l-ph 3 T 2-ph Rot-ph 

= A e " (12 K)/T + BT 7 + C T 3/2 e _(8 ' 6 K)/T , (1) 

where the first term comes from one-phonon absorption, 
with A = 130 s -1 [55]; the second from 2-phonon scat¬ 
tering (one phonon absorbed and another emitted), with 
B = ( 8.8 and 7.6) x 1CT 3 s " 1 K " 7 for 0.6 K and 1.0 K, 
respectively; and the third from roton-phonon scattering 
(a roton absorbed followed by a phonon emitted), with 
C = 18 s - 1 K -3 / 2 . The A, B and C coefficients were 
calculated for T < 1 K, however, their temperature de¬ 
pendences are weak compared to the overall sizes of the 
terms. In Fig. |T[ a plot of these terms up to 2.2 K is 
shown to illustrate their expected size. 

Experiments that have studied the temperature- 
dependence of r^ p 1 (T) have done so up to 1.2 K [351 1531 
155] , 1.5 K [32], and 1.6 K JBU]. In all these studies, us¬ 
ing r~p(T) = t^_ p h only was sufficient to describe the 
data. In the second part of this paper, we study UCN 
up-scattering to 2.4 K. A detailed data-fitting procedure 
is used to test the validity of the expected temperature- 
dependence of T~p(T). 


II. THE EXPERIMENTAL SETUP 

The general layout of the experiment installed on the 
PF1B CN beam line [5T] in 2011, the details of which 
are described later, is shown in Fig. [2] A continuous 
flow 3 He refrigerator with a measured cooling power of 
60 mW at 0.6 K [52] is situated in the main cooling 
tower. This refrigerator is thermally connected, via a 
copper heat exchanger, to a horizontal stainless steel tube 
filled with superfluid helium with 50 mm diameter and 
45 cm length. This filled-tube supplies superfluid helium 
to the converter vessel and also serves for heat trans¬ 
fer. UCNs produced in the helium by the CN beam are 
extracted vertically from the back of the converter vessel 
and guided by tubes in a “fl”-shaped geometry to a UCN 
detector at room temperature. 

The superfluid 4 He used for the converter is purified 
through a superleak made from AI 2 O 3 powder with a 
grain size of 50 ran. The fabrication technique is de¬ 
scribed in Ref. (53]. Superfluid flow through the su¬ 
perleak is induced via the thermomechanical pumping 
(or fountain effect) with a heater, similar to that of 
Ref. [53]. The 3 He to 4 He isotopic purity from this pro¬ 
cess is < 4 x 10 -10 [55]. 

The incoming CN beam first passes through front alu- 



FIG. 1. (Color online) The sizes of the three main processes 
contributing to the up-scattering of UCNs by excitations in 
superfluid helium from Eq. 0 . derived in Ref. [46]. The thick¬ 
ness of the 7 © p h line covers the range of values for T < IK 
due to the different values of B (see text). The dotted lines 
are used to indicate that the up-scattering rates are only ap¬ 
proximate for T > 1 K due to temperature-dependences in 
the A, B and C coefficients. 



FIG. 2. (Color online) The experimental installation at the 
PF1B CN beam line. For clarity, extraneous elements of 
the original photograph are erased. The added labels indi¬ 
cate: © main cooling tower; © horizontal helium connection; 
@ converter volume; © mechanical UCN flap valve actuator; 
© UCN extraction guides; and © UCN detector. The direc¬ 
tion of the CN beam is shown as the blue arrow. Additional 
lead shielding surrounding the apparatus was added later. 







3 


minum windows of the vacuum vessel, the two radia¬ 
tion shields, and the superfluid containment vessel be¬ 
fore reaching the converter through a 1 mm thick beryl¬ 
lium sheet, which acts as the front window for CNs and 
the front wall of the converter vessel for UCNs (Fig. [3]). 
The side walls of the converter vessel are made from four 
beryllium-coated aluminum plates assembled together 
with aluminum screws (Fig. [4]). This assembly is clamped 
onto a 8 cm long beryllium-coated copper end-section, 
where the mechanical UCN valve is located. Another 
1 mm thick beryllium sheet serves as the back CN win¬ 
dow and UCN wall. The internal cross-sections of the 
converter vessel are 8x8 cm 2 in the section made of the 
aluminum plates and 10 x 10 cm 2 in the copper end sec¬ 
tion. The volume of the converter vessel with the UCN 
valve closed is Konv = 4.0 liters. 



FIG. 3. (Color online) Cross-section of the superfluid helium 
converter volume and UCN extraction system. © Isolation 
vacuum window; © outer radiation shield window; © in¬ 
ner radiation shield window; © superfluid containment ves¬ 
sel (shown transparent); © front beryllium sheet; © beryl¬ 
lium coated aluminum side walls; © beryllium coated copper 
end section; © beryllium UCN flap valve; © back beryllium 
sheet; © pulley of UCN flap valve mechanism; © polished 
stainless steel conical guide; © polished stainless steel UCN 
guide; © thermal anchoring to 2nd stage radiation shield; and 
© UCN guide adapter (50 mm to 66mm ID). The directions 
of the incoming CN beam (blue arrow) and the remaining 
UCN guides to the detector (red arrow) are also shown. (The 
concentric cylindrical side-walls of the radiation shields and 
vacuum vessel are not shown.) 


The beryllium coatings on the aluminum plates and 
copper were produced at the Petersburg Nuclear Physics 
Institute, Russia, by performing two sputtering proce¬ 
dures with intermediate cleaning to reduce defects of 
coating caused by dust on the surface. Each is estimated 
to be approximately 400 nm thick. 

Measurements of the converter temperature are made 
using 3 calibrated Cernox™ sensors from Lake Shore 
Cryotronics placed outside the converter volume but in¬ 
side the superfluid containment vessel: one is placed at 
the front of the converter vessel and two at the back on 



FIG. 4. (Color online) View of the converter vessel before the 
superfluid containment vessel is mounted. The added labels 
and arrows indicate: © screws for assembling the side walls; 
© claw clamps for mounting the front sheet to the side walls; 
© claw clamps for mounting the side walls to the copper 
end-section; and © wiring for thermometer at the front end 
of converter volume. 


the copper end-section at different heights. During the 
measurements, the temperatures read at these three sen¬ 
sors were the same to within < 20 mK. 

The UCN flap valve (Fig. [5]) consists of a circular 3 mm 
thick beryllium disc connected to an axle of a pulley sys¬ 
tem. A stainless steel cable, driven by a pneumatic piston 
located outside of the isolation vacuum, is used for clos¬ 
ing the valve. A spring, and the weight of the plate itself, 
provides the force for opening it. The stainless steel axle 
and wire supporting the plate inside the converter ves¬ 
sel were coated with nickel (Uni = 252 neV) via brush 
electroplating. 



FIG. 5. (Color online) The mechanical UCN flap valve viewed 
from the incoming direction of the CN beam. The axle and 
wires supporting the flap can be seen, as well as the exit 
for UCNs at the bottom of the stainless steel conical guide 
section. 
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The first section of the vertical UCN extraction guide 
system is the polished stainless steel ( Uss = 185 neV) 
conical guide piece. It reduces the radiative heat transfer 
into the converter volume from the guides and also aids in 
aligning the outgoing UCN trajectories vertically so that 
they can overcome the gravitational potential with fewer 
reflections off the walls of the guides. The diameter of 
the entrance to the cone is 23 mm. The rest of the UCN 
extraction guides, also made of polished stainless steel, 
has a 50 mm ID until around half-way on the horizontal 
section after which it increases to 66 mm. The vertical 
distance from the closed UCN valve to the bottom of the 
horizontal section is 28 cm. 

The horizontal guides have a total length of « 80 cm 
and the vertical drop to the UCN detector is 1.2 m. The 
detector is a 3 He gaseous wire chamber with an aluminum 
(Uai = 54 neV) front window. The dead time of this 
detector was measured to be 2 /is [4ll . 

The conductive heat load down the walls of the guides 
is reduced by thinning out the straight sections of the 
guides to a wall thickness of 0.5 mm. In addition, the 
guides are thermally anchored to the inner and outer 
radiation shields at some points. A difficulty of using 
the window-less extraction scheme is the heating due to 
thermal radiation. Radiation baffles can not be used 
since they would obstruct UCNs. To reduce this radia¬ 
tive load, the thermal anchoring to the inner heat shield 
is done both at the copper block and at the 90° bend 
just above it to reduce the temperature of surfaces that 
have a direct view to the converter. In the current setup, 
there is no mechanism for suppressing superfluid helium 
film flow. With a similar extraction guide system and a 
3 He refrigerator with similar cooling power, a converter 
temperature of < 0.72 K has previously been achieved 

mmm- 

The PF1B CN beam line is a 74 m long m = 2 super¬ 
mirror guide coupled to a liquid deuterium cold source 
[SEJ- The cross-section of the neutron guide in the case¬ 
mate has a dimension of 6(W) x 20(H) cm 2 . This was 
cut with an aperture so that two extension guides with 
cross-sections of 12 x 5 cm 2 , forming a combined length of 
2 x 2 m, could be used. The beam was then further cut to 
8x5 cm 2 just before the converter vessel. The total vol¬ 
ume in contact with the beam (neglecting its divergence 
within the converter vessel) is Vqn = (2,300 ± 100) cm 3 . 
A CN beam monitor detector behind a small circular 
aperture was placed behind the converter volume outside 
the cryostat for monitoring the relative CN flux. 

An older characterization of the PF1B beam inside 
the casemate m measured a total capture flux of 1.4 x 
10 10 cm -2 s _1 , with a spectrum that peaks at 4 A with 
a differential capture flux of 2 x 10 9 cm -2 s _1 A , and 
at 8.9 A, where single-phonon UCN production takes 
place, a differential capture flux of 7 x 10 8 cm -2 s _1 A 
In around 2007, after the measurements of Ref. [61] . 
parts of the guide were repaired and the in-pile part 
upgraded, which increased the total capture flux to ~ 
2.2 x 10 10 cm -2 s _1 [67]. The capture flux at 8.9 A may 


have changed more strongly due to colder neutrons mak¬ 
ing on average more reflections on the guides. The exten¬ 
sion guides, despite accepting full divergence of the pri¬ 
mary beam, would have also introduced some loss. The 
estimated flux for UCN production in the converter is 
only expected to be within ~ 30% of the above values 
m ■ Since this is used only for comparing the measured 
UCN production rate in this paper, it is considered to be 
of sufficient precision. 

While the mean-free-path of 8.9 A neutrons in super- 
fluid helium is ~ 17 m, it decreases to ~ 70 cm at 3 A [68], 
the center of the broad multi-phonon production peak. A 
two-component epoxy doped with boron (natural isotopic 
composition) carbide powder was coated on the cylindri¬ 
cal side-walls of the aluminum radiation shield surround¬ 
ing the converter volume to absorb scattered neutrons. A 
lead wall at least 15 cm thick in all places was sufficient 
for biological shielding. 

The background rate at the UCN detector when the 
CN beam is switched on comes primarily from prompt, 
scattered CNs. An upper-limit, deduced from having the 
converter at 2.4 K and UCN valve closed, is < 12 s _1 . 
This rate is higher than typical backgrounds when using a 
monochromated 8.9 A CN beam, common for superfluid 
helium UCN production. The ambient background rate 
when the CN beam switched off is < 0.05 s _1 . Further 
details of the experimental setup and the data described 
in the subsequent sections can be found in Ref. [69 ]. 

III. THE UCN POPULATION DURING 
MEASUREMENTS 

An accumulation measurement starts with the UCN 
valve closed and no UCNs in the converter volume. The 
CN beam is switched on at t = 0. Assuming the UCNs 
are in mechanical equilibrium [70], if N(E, T, t) d E is the 
number of UCNs inside the volume with total energy 
between E and E + d E, the build-up of UCNs will follow 

N(E,T,t) = N oa (E,T) {1 — exp[—f/Ttot(-E, T)]} (2) 

where 

N x (E,T) = P(E, T ) U C n not (E,T). (3) 

P(E,T) is the differential UCN production rate per unit 
time and volume. 

The total characteristic storage time constant, 
r tot (E, T), can be written in terms of loss rate constants 

r~l (E, T) = Tq 1 (E) + r" 1 (T) + + r^ 1 , (4) 

where t^ 1 (E) consists of: (1) capture or up-scattering 
losses through interactions with both the vessel walls 
(due to the UCN wave-function tunneling into the bulk 
material of the walls) and impurities frozen on the walls, 
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and (2) UCNs escaping through gaps between walls of the 
vessel and through the imperfect seal of the UCN valve 
when it is closed. Tq(E) is higher for larger E due to an 
increased wall collision rate, as well as an increase in the 
loss probability per reflection 123 ■ The loss due to up- 
scattering of UCNs by excitations in superfluid helium, 
r^ p 1 (T), is described in Eq. [lj The loss of UCNs on im¬ 
purities distributed in the converter, which is primarily 
due to remaining 3 He contamination after super-leak fil¬ 
tration, is denoted by t~^ s and is E- and T-independent. 
Finally, Tp is the neutron /3-decay lifetime. 

P(E,T) is dependent on the converter temperature T 
(see the appendix). The super-thermal down-conversion 
process is expected j47j to fill phase-space with a con¬ 
stant density so that P(E, T) <x \[E. It should be noted 
that, since T tot (E,T) is shorter for higher UCN energies, 
the spectrum of accumulated UCNs N oa (E, T) is not the 
same as the one produced via P(E,T). 

After reaching saturation in the accumulation, i.e., 
t T to t{E,T) for all E of UCNs storable in the con¬ 
verter volume, the next step is to switch off the CN beam 
and then open the UCN valve. Denoting the time when 
this occurs as t = t 0 , then the open vessel will empty 
according to 

N(E, T, t ) = N(E, T, t 0 ) exp[—(t - t 0 )/r e (E, T)} , (5) 
where 

t~\E,T) « r t ~ 0 l(E,T) + t-\E) (6) 

is the emptying rate constant. t^(E) is the rate con¬ 
stant of UCNs escaping through the exit of the opened 
valve. This is only approximate since it assumes that 
opening the UCN valve does not change t 0 (E). Also, 
since r e (E,T) is shorter for higher E, the spectrum 
N(E,T,t ) evolves with time. 

For a continuous measurement, which also starts with 
no UCNs in the converter volume, the CN beam is 
switched on with the UCN valve left open. Then the 
build-up of UCNs (the term accumulation is reserved 
for when the UCN valve is closed) will follow Eq. [2] but 
with r e (E,T) in place of Ttot(E,T). After saturation is 
reached, i.e. t r e (E,T), the CN beam is switched off 
(at t = to)- The population of UCNs in the vessel will 
empty according to Eq. [5] as before. 


IV. 1.08 K ACCUMULATION MEASUREMENT 

An accumulation measurement with the converter tem¬ 
perature at T = 1.08 K was performed. The count rate 
observed at the UCN detector during this measurement 
is shown in Fig. [6] Due to poor cooling of the radia¬ 
tion shields surrounding the volume this was the lowest 
temperature attained in the experiment. 

Due to the UCN valve not being fully tight against 
leakage of UCNs when closed, UCNs can reach the detec- 



FIG. 6. (Color online) UCN count rate during the 1.08 K 
accumulation measurement. At 187 s (= to) the UCN valve 
was opened and CN beam switched off. At 258 s the UCN 
valve was closed again, before all the UCNs were emptied from 
the converter volume. The fits of the build-up and emptying 
curves, details of which are described in the text, are shown 
as solid lines. 


tor during the accumulation period. The detected rate, 
denoted by N acc (T, t), is related to the number of UCNs 
in the converter vessel N ( E , T, t) from Eq. [2j Since the 
UCN detector does not provide energy information, it 
only observes the integral over E for UCNs storable in 
the vessel. Thus, N acc (T,t) is given by (neglecting the 
time delay due to the time-of-flight of the UCNs) 

N acc {T,t) = J e m (E)T{^(E)N 0O (E,T) 

x {1 - exp[—t/r tot (U, T)]} d E , (7) 

where e ex (E) is the efficiency of the UCN extraction and 
detection, and 'r k / a l k (U) is the rate constant of UCNs es¬ 
caping through the leak at the UCN valve. It will be seen 
later that r^ k (E) -C t^I(E), so the size of leak does not 
affect Tt 0 t(E,T). N acc (T,t) can be used for monitoring 
the relative number of UCNs inside the vessel. TE~ k (E) 
and the absolute value of N acc (T, t) are not important 
in this analysis. The limits of integration over E are 
from the minimum UCN energy that can reach the de¬ 
tector due to vertical extraction to the maximum energy 
of UCNs storable in the converter and transportable by 
the guides. 

When the UCN valve is opened and the CN beam 
switched off, the accumulated UCNs start emptying from 
the volume. Using Eq. [5] (similar to Eq. [7]), the count rate 
at the detector during emptying (t > to) is given by: 

iVe mp ty(T,i) = J e ex (E)T~ 1 (E)N(E,T,t 0 ) 

x exp[-(t - t 0 )/r e (E,T)] dE . (8) 
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In order to fit the data without a priori knowledge 
of the energy dependence of 7 tot (E,T) and r e (E, T), the 
integrals of Eqs. [7] and [ 8 ] can be approximated by a dis¬ 
crete sum of terms with different amplitudes a.i and fa 
and time constants T to t,i and r e ^ given by 

V acc 

N acc (T,t ) ~ ^ Ui {1 - expl-t/Ttot.iCO]} (9) 

i =1 

and 

'H empty 

N empt y{T,t) SS ^2 fa exp[—(t — to)/Te t i(T)] . (10) 

i—1 


The size of rj acc or ?7 em pty required to satisfactorily de¬ 
scribe the data (i.e. xt ~ 1) depends on the strength of 
the E-dependence in r tot (E,T) or r e (E, T), as well as on 
the statistical precision of the data. 

Using rja.cc = 1 is sufficient for the fit of iV acc (1.08 K, t) 
from 0 < t < t 0 using Eq. [9j which gives xt = 0-80 with 
r to t(1.08 K) = (20.3 ± 1.2) s. Due to rj acc = 1 being a 
valid approximation, we drop the i-subscript and write 
simply T t ot. 

Using 77 empty = 2 to fit N em pty (1.08 K, t) from t > t 0 
to when the UCN valve was closed gives xt = 11, with 
the values: fa = (5,060 ± 160) s _1 , r e ,i = (6.1 ± 0.2) s, 
fa = (2,530 ± 180) s _1 , and r 6i2 = (16.3 ± 0.4) s. For 
comparison, using 77 em pty = 1 gives xt = 20. The need 
for ?7empty > Vacc is due to the increased statistics of the 
emptying data and the larger E-dependence of r e (E,T) 
(see Eq. [6]). 

Following previous papers 53 ESj, we write Eq. [3] in 
its integral form since it is ensemble averages we are con¬ 
cerned with. The total number of UCNs accumulated in 
the vessel at saturation can thus be expressed as 


is the probability a UCN can find the valve exit before 
becoming lost in the vessel. 

The discrete ?7 em pty = 2 approximation of the empty¬ 
ing curve (Eq. 10) can be applied to Eq. 12 to give 


Ncmpty (-C) ^ex,l € v ,l IVoo,! + £ex,2 Cv,2 l^oo,2 , (14) 


where 6 v,z — [^tot "re,i]/r"tot, ^oo.i — 
and Noo = AA^i ± A/*,^. If we define an effective extrac¬ 
tion and detection probability value 


^ex 


^ex, 1 ^v, 1 iVoo,l + ^ex ,2 ^v ,2 A/qo ,2 

£v,l Noc,l ± £v ,2 


(15) 


then we can write 


A/e mpty{T) — Cex^v,! N<x,l + £ v ,2 ^ 00 , 2 ) ■ (16) 


e ex is the weighted average of e ex ,i and e e x ,2 with weights 
£v,i N<x>,i and e v . 2 Noc, 2 , respectively. Similarly, we can 
further define 


e v = 


Cv,] A^oq^i T e v ,2 JVoo.2 
IVoo,l + iVoo .2 


(17) 


so that 


A^empty (T) — £ e x ^v{^oo,l + ^ 00 , 2 ) ■ (18) 

From combining Eqs. [I 8 | and [12} the volumetric UCN 
production rate, without making corrections for loss 
mechanisms, can be deduced to be 

e ex e v (1.08 K) E(1.08 K) = (1.96T0.14) cm " 3 s ’ 1 . (19) 


If we wish to correct for the loss due to e v (E,T), we 
put the values for e v i used in Eq. 14 into Eq.[l 6 |and then 
combine with Eq. [12] to get 


N 00 (T) = P(T)V CN T t0t (T). (11) 


Cex P(1.08 IC) = (6.9 ± 1.7) cm -3 s " 1 . (20) 


This expression can be used to indirectly determine the 
volumetric UCN production rate P(T). 

Summing the counts observed in the emptying curve 
gives A/empty (1-08 K) = 91,700 ± 300. To extract P(T), 
A^empty (T) is used to estimate the number of UCNs 
present inside the converter volume at the end of ac¬ 
cumulation, N(E,T,to). If we assume that saturated 
accumulation was reached, i.e. N(E,T,t 0 ) = Nqc^EjT) 
for all E, then the two quantities are related by 
/»00 

Afempty CO = / N empty (T,t) dt ( 12 ) 

Jt 0 

nOO 

= / e e x{E) e v (E, T) N ao {E, T) dE, 

Jto 

where 


= r e (E,T) _ r tot (E,r)-r e (E,r) 

_ t v (E) T to t(E,T) 


( 13 ) 


The correction for e e x is more difficult to make since it 
can not be derived from observed quantities in this ex¬ 
periment, and it depends on the UCN spectrum. UCN 
transport and detection is however not expected to be 
strongly E-dependent due to the few reflections (com¬ 
pared to storage) and the vertical drop before the de¬ 
tector. The detection efficiency and transmission at the 
window our UCN detector is ~ 80%. Typical transmis¬ 
sion per unit length of polished stainless steel UCN guides 
of this diameter are ~ 90 %/m 70j. And in a previous 
experiment with a similar UCN extraction guide geome¬ 
try and UCN detector of the same type, from comparing 
with in situ CN beam characterization measurements, a 
value of e ex ~ 38% was obtained [55] . 

The PF1B volumetric UCN production rate was calcu¬ 
lated to be (13.9 ±0.9) cm -3 s ” 1 in Ref. [7T] by using the 
beam characterization data from Ref. |61j and an upper 
energy cut-off of U^, e — U^ e , where En e = 18.5 neV is 
the optical potential of the converter. This value has no 
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adjustment for vertical extraction. 

When UCNs leave the helium, they get a boost in ki¬ 
netic energy of U Be . The height of the vertical extraction, 
measured from the flap valve to the horizontal guide, is 
h ex = 28 cm, which corresponds to a potential difference 
of U ex = 29 neV. Assuming a constant phase space den¬ 
sity production spectrum, i.e. P(E, T) oc VE, the correc¬ 
tion factor of the production rate from Ref. 1711 for verti¬ 
cal extraction is: [(U Be - U Be ) 3/2 - (U ex - U He ) 3/2 ] /(U Be - 
U Be ) 3 ^ 2 = 0.96. This corresponds to a volumetric pro¬ 
duction rate of P = (13.3 ± 0.9) cm -3 s _1 . 

In our setup the UCN extraction guides are stainless 
steel so it is not expected that all UCNs with E up to 
U Be outside the converter make it to the detector. (UCNs 
with E > Uss can only do so if they only reflect off the 
UCN guides above the level of the converter, or if they 
make only glancing reflections off the guides only.) A 
pessimistic case is that all UCNs with E > Uss outside 
the converter are lost, then the correction factor becomes: 
[(U SS - U He ) 3/2 - (U ex - C/ H e) 3/2 ]/(UBe - U He ) 3/2 = 0.67, 
corresponding to P = (9.3±0.6) cm _3 s _1 . These values 
do not take into account the variation of the actual CN 
flux compared to that from Ref. m as mentioned earlier. 



FIG. 7. (Color online) Continuous measurements for selected 
temperatures to demonstrate the changes in the build-up, sat¬ 
uration, and emptying. The duration the CN beam is left on 
is not the same between the measurements. The last 50- 
60 s before the CN beam is switched off is used to determine 
N C (T). 


V. 1.2 K TO 2.4 K CONTINUOUS 
MEASUREMENTS 


Continuous measurements over a range of converter 
temperatures were used to study r up 1 (T). The UCN 
count rates for selected measurements are shown in 
Fig. 0 and the saturated count rates N C (T) as a function 
of temperature are plotted in Fig. [8] The temperature 
during each measurement was stable to < 10 rnK. 

N C (T) is given by the rate of continuously produced 
UCNs in the converter vessel that find the exit of the 
valve (Eq. [l3]) and get detected; i.e., 

N C (T) = f P(E , T) Vcn e ex (E) ^ [ ^ d E . (21) 

J T e (E,l ) 


Since all the temperature-dependence is isolated in 
r^ p 1 (T), an expression appropriate for fitting N C (T) to 
study its temperature-dependence can be written as 


N C (T) = 


P(T)V CNCe 


r-1 


N C (T ) 

>rel(T) 


T 0 + Tu P (T) + T abs + Tp + Tv 

Po Ucn CexAT 1 


-1 


r" 1 -hr" 1 

r e,T->0 ^ Tu P 


(T)’ 


( 22 ) 


which comes from expressing the integral of Eq. ^T] with 
the ensemble averaged values e ex , A7 1 , and Tq" . We 
do not require to know these values, they need only to 
be constant with temperature. All the rate constants 
contributing to t~ 1 besides rJ p 1 (T) are combined into 
T eT^ 0 , which is used as one of the two free parameters of 
the data fitting. The other free parameter is the scaling 






FIG. 8. (Color online) Bottom : The saturated count rates 
from continuous measurements N c between 1.16 - 2.40 K. The 
region around T\ is shown zoomed-in. The solid green line 
is the fit for T < T A using r~p(T) = T^ h . N C (T)/P re i(T) 
is being fitted in the analysis, but N C (T) is shown here. See 
text for details. Top : The relative error of each N C (T)/P le i(T) 
point and the residuals (fit minus data) of the r^p (T) = ph 
fit are also shown. T\ is marked by the vertical dashed lines. 













































factor PoVcNCexT” 1 , where Pq = P(T > 0 K). The 
T-dependence of production is contained in the relative 
UCN production rate P re i(T), i.e. P re i(T) = P(T)/P 0 . 

P re i(P), required for the data-fitting, is calculated via 
the van Hove relation m using experimental data from 
Ref. [73] as described in the Appendix. The results are 
shown in Fig. [9j P re i(T) decreases with temperature as 
T\ = 2.172 K is approached and then increases. Qualita¬ 
tively, the kink in the N C (T) data at T\ can be explained 
by Prei(P)- Accurately describing N C (T > T\) is however 
difficult due to effects such as helium vapor bubbling in 
the liquid, and is not attempted here. 


routine was varied to ensure we are finding the global xl 
minimum. 

When T < T\ is used the best fit is obtained using 
r~p(T) = t^_ p h producing xl = 3.8. A plot of this fit 
and its residual is shown in Fig. [ 8 ] The next best fit 
comes from using T“ p 1 (T) = r” 3 h + r”^ h with xl = 8 T. 
The former fit, according to the xl distribution (z/ = 30) 
m is 10 23 times more likely to be a better model for 
the data. Any fit that includes T rot '_ ph , which should be 
comparable in size at this temperature range, is much 
poorer. 



FIG. 9. (Color online) The relative UCN production rate 
Prei(T) calculated using inelastic neutron scattering data and 
an analytic approximation of the PF1B CN beam spectrum 
(see the Appendix for details of the calculations). The solid 
green line is the interpolation between the points with Eq. [23] 
used in the data analysis. 


For the data analysis, interpolation between the 
P re i(T) values is required. Below T\. with no known 
physical significance, the expression 

P re i(T < T X ) {1 - exp [{T-p 2 )/p 3 }} (23) 

is used to fit the P re i(T) values for this purpose. This 
gives pi = 0.999 ± 0.003, p -2 = (2.74 ± 0.10) K, and 
p 3 = (0.22 ± 0.04) K, with a xl = 0-62. 

The uncertainty in the N C (T) / P re \{T) points used for 
fitting with Eq. [22] comes from combining in quadrature 
the statistical uncertainty in N C (T) (0.2-0.8%) and the 
uncertainty from the interpolation of P re i(P) (0.3-3%). 
This is shown in Fig. [H] 

To determine the best form of that describes the 
data, non-linear fitting with different combinations of the 
terms for T~p{T) from Eq. ^is done. A comparison be¬ 
tween fitting with the fixed coefficients A = 130 s -1 , 
B = 7.6 x lO” 3 s _1 K" 7 , and C = 18 s" 1 K" 3 / 2 for 

T 2 -phi an d T rot-ph’ respectively, is shown in Table I In 
all the analyses, the starting values of the optimization 


TABLE I. Fitting N C (T < T\) and N C (T < 1-95 K) with 


Eq. 22 for different combinations of r up and using fixed co¬ 


efficients A, B, and C according to theoretical predictions 
(see discussion around Eq. [l]). The best-fit parameters and 
reduced chi-squared xl = 30 and 19 respectively) are tab¬ 
ulated. Bold font indicates the best fit. 


7- 1 

up 

FbVcNCexTV 1 
[s- 2 ] 

r e,r-K> [ S ] 

xl 

T < T\ (y = 30) 

T -1 
* 1-ph 

245 ±1 

26.25 ±0.10 

31 

T 2-ph 

680 ±2 

10.56 ±0.05 

3.8 

rot-ph 

547 ±1 

12.77 ± 0.05 

56 

— 1 —1 

^2-ph ' ^rot-ph 

1229 ±3 

5.77 ± 0.02 

18 

— 1 , —1 

T l-ph + T 2-ph 

927 ±2 

7.52 ±0.03 

8.1 

-1 , -1 

T l-ph ' 'rot-ph 

792 ±2 

8.59 ±0.04 

46 

— 1 , — i . -1 

T l-ph T 2-ph T rot-ph 

1475 ±4 

4.73 ±0.02 

19 

T < 1.95 K (u = 19) 

T i- P h 

251 ±1 

25.36 ±0.11 

11 

^2-ph 

686 ±2 

10.44 ±0.05 

1.7 

r -l , 

rot-ph 

564 ±1 

12.19 ±0.05 

25 

— 1 , -1 

2-ph ' ^"rot-ph 

1250 ±3 

5.62 ±0.03 

8.0 

-1 . -1 

T l-ph + T 2-ph 

937 ±2 

7.39 ±0.03 

3.0 

- 1 , -1 

T l-ph ' T rot-ph 

815 ±2 

8.24 ±0.04 

18 

— 1 . — i , -1 

T l-ph T 2-ph T rot-ph 

1501 ±4 

4.60 ±0.02 

7.7 


As can be seen from the fit residual shown in Fig. [ 8 j 
a large contribution to the xl comes from T > 1.95 K. 
This is true for the other fits as well. At these higher 
temperatures, approximations used in the fit break down 
on the ~ 10% level. This indicates there are new phe¬ 
nomena causing either an increase in or r^ -1 , or 

reducing P re i(T) beyond that predicted by the calcula¬ 
tions with neutron scattering data. Since these new ef¬ 
fects are not included in the description provided thus far, 
it is more valid to restrict the data fitting to T < 1.95 K 
only. The results of repeating the previous analysis in 
this temperature range is also shown in Table |T[ The 
same conclusion is reached: the best fit is still from 
r^p 1 = t^_ p h , with the xl decreasing to 1.7, while the 
rest of the xl values decrease by a factor of ~ 3. The 
probability that ph is a more likely fit than 

r~p{T) = Tfph + t^_ p h is now 3000 times greater only 
m- In all subsequent analyses, we restrict to using 
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T < 1.95 K. 

Also explored are potentially incorrect fixed values for 
A , B, and C affecting the quality of the fits. To do this, 
the values were varied to see their impact on the y 2 . 
This is done by allowing these values to be optimized 
by the Levenberg-Marquardt non-linear regression data- 
fitting algorithm used in the analysis so far. However, if 
the y 2 dependence of a certain variable is weak, then the 
algorithm does not converge. In this case, the variable 
was varied manually. 

First, we make the fit using r^ p 1 (T) = r^ ph and vary 
the parameter B. This had to be done manually since 
the y 2 does not change for a wide-range of B values, 
as shown in Table [nj y 2 stays the same because the 
values of t b t^o an d -Fb^CNCexT ” 1 are able to compen¬ 
sate, however, at some stage r ei T->-o becomes unphysi¬ 
cal when compared to the observed 7 +, (1.08 K) values 
from the accumulation measurement of ( 6.1 ± 0 . 2 ) s and 
(16.3 ± 0.4) s (see previous section). These values are 
expected to be similar to those when T — > 0. We thus 
make the assertion that r e> T-y o should fall in the range 
conservative range of (5 - 20) s. This leads to a value for 
B = (4 - 16) x 10 ~ 3 s _1 K -7 being sufficient to describe 
the data. The fitted PoVcNEexAT 1 values can be com¬ 
pared to the expected value too, however, this is much 
less stringent due to the number of parameters required 
to be estimated. 

TABLE II. Fitting N C (T < 1.95 K) with Eq. || using 
T up {T) = t^_ p h and with B varied manually. Bold font in¬ 
dicates fits that satisfy the acceptable limits for r e ,T ^o (see 
text). 


B 

[icr 3 s- 1 K- 7 ] 

0.76 

4 

16 

76 

xl 

1.7 

1.7 

1.7 

1.7 

EbVcN e exTV 

[s- 2 ] 

69 ±0.2 

361 ± 1 

1443 ± 4 

6860 ±20 

T e ,r-K) [s] 

104 ±0.5 

19.8 ±0.1 

4.95 ±0.02 

1.044 ±0.005 


Next we analyze the effect of varied coefficients when 
Tup 1 CO = T 2 -lh + r rot- P h is use d for the fit. Like pre¬ 
viously, B has to be varied manually, however C can 
be optimized with the fitting algorithm. The results of 
this analysis are shown in Table m In the optimiza¬ 
tion, the values of C are restricted to C > 0. It can be 
seen that the optimized C values for all values of B used 
are consistent with zero. Because of this, the optimized 
T e ,T—>-o and y 2 are similar to those from the fits with 
r~p(T) = p h and B fixed. The weak dependence of 
the fit on B, and the previously quoted acceptable range 
for B 1 are again confirmed. From this analysis, we de¬ 
duce the largest possible value for C is < 5.3 s -1 K -3 / 2 
(95% C.I.) in order to keep r e> T->o ~ (5 - 20) s, which 
equates to B = (4 - 16) x 10 -3 s _1 K -7 , the same range 
from the previous analysis. 

Finally, we vary all three coefficients in the r^ p 1 (T) = 


Wph + r 2 ^ph + T r “ t _ ph fit. The value of A also has to be 
varied manually; as before, B has to be varied manually; 
and C can be optimized algorithmically. The results of 
this analysis are shown in Table |IV| Both A and B have 
been varied by two orders of magnitude or more but all 
the optimized C values are still consistent with zero. For 
all values of B 1 the smaller the A value the lower the y 2 . 
This further reinforces the absence of r r “ t . ph in the data, 
as well as demonstrates the absence of Tj 7 ph . 

Since this analysis is the joint variation of two param¬ 
eters (A and B) with the other 3 parameters (r e j_, 0 > 
-PoVcNfexAT 1 ) and C) optimized by the Levenberg- 
Marquardt algorithm, we have separated the y 2 fitting 
problem into two parts: a variation of y 2 for the fit of 
the data to 3 parameters (with v = 18) and a varia¬ 
tion of Ay 2 about the minimum y 2 with 2 parameters 
(v = 2). These two variations separately follow their 
appropriate y 2 -distributions |74j. The limits of 95% 
C.I. are thus given by contours of min(y 2 ) + 6.14 (or 
min(y 2 ) + 0.34 = 2.17). 

The limit of A is found by finding the widest part 
of the 95% C.I. contour that also satisfies the restric¬ 
tion that o = (5 - 20) s, which again equates to 
B = (4-16) x 10 -3 s - 1 I<+'. This occurs for B = 
14 x 10 ~ 3 s _1 K -7 , and sets the limit A < 61 s _1 (95% 
C.I.). The limit on C requires having A optimized for 
each fit which occurs for A —> 0. This provides the 
limit of C < 5.3 s _1 K -3 / 2 (95% C.I.) from this analysis, 
the same as the limit from the previous analysis. 

VI. CONCLUSION 

This paper presented UCN production measurements 
using a superfluid 4 He UCN source with a 4-liter 
beryllium-coated vessel, mechanical UCN accumulation 
valve, and vertical, window-less stainless steel extraction 
system. Experiments were performed at the PF1B CN 
beam at the ILL. 

In the first part, with an accumulation measure¬ 
ment at 1.08 K an uncorrected UCN production rate 
of e ex e v (1.08 K)P(1.08 K) = (1.96 + 0.14) ernes ’ 1 
was observed. After correcting for e v (1.08 K), whose 
value is dominated by a poor r to t(T), a production rate 
of e ex P(1.08 K) = (6.9 ± 1.7) cm _ 3 s _1 was extracted. 
The approximate size of e ex estimated from other ex¬ 
periments is ~ 40%. The calculated production rate is 
P ~ (9-13) cm -3 s _1 , where the range comes from de¬ 
tails in the UCN extraction process. Furthermore, these 
rates are based on older beam characterization measure¬ 
ments EH and could be different to the actual flux in the 
experiment by ~ 30%. 

The short r to t(l-08 K) = (20.3+1.2) s might be caused 
by a combined effect of UCN losses due to imperfec¬ 
tions of the beryllium coatings of the walls and patches 
with high losses on the complex UCN valve mechanism, 
which both can be improved. Better cooling of the con¬ 
verter will improve things further, since T up (1.08 K) « 
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TABLE III. Fitting N C (T < 1.95 K) using Eq. 22 and t up 1 (T) = t 2 _ ph + i" ro /_ ph with B varied manually and C optimized 
algorithmically. Bold font indicates solutions that satisfy the acceptable limits for r e ,T->o (see text). 


B [10 -3 s —1 K —7 ] 

0.76 

4 

7.6 

16 

76 

xt 

1.8 

1.8 

1.8 

1.8 

1.8 

C [s —1 K —3 / 2 ] 

(1 x 10" 4 )±0.1 

(1 x 10 ~ 3 ) ± 0.7 

(3 x 10 " 4 ) ± 1.3 

(4 x 10 ~ 3 ) ± 2.7 

(2 x 10" 2 )± 1.3 

T'e.T-t0 [S] 

104 ±6 

20 ±1 

10.4 ± 0.6 

5.0 ± 0.3 

1.04 ±0.06 


TABLE IV. Fitting N, (T < 1.95 K) with Eq. 22 using r^, 1 (T) = r^ h + r^ ph + T~/_ ph with A varied manually, B varied 
manually, and C optimized algorithmically. Bold lont indicates solutions that both satisfy the acceptable limits for r e ,T->o as 
well as the 95% C.I. limit of xt (see text). 




1.3 

18 

33 

61 

130 

1 

W 

1 

cn 

CO 

1 

O 

0.76 

xt 

C [s^ 1 K" 3 / 2 ] 

As,T—M) [s] 

2.0 

(1 x 10 -4 ) ± 0.1 
100 ±6 

3.7 

(3 x 10~ 5 ) ±0.2 
66 ±4 

4.8 

(6 x 10" 5 ) ± 1.4 
51 ± 3 

6.3 

(5 x nr 5 ) ± o.4 
36 ±2 

8.1 

(1 x 10 -5 ) ± 0.8 
20 ±1 

4 

xt 

c [s- 1 k- 3 / 2 ] 

Te,T -¥0 [s] 

1.9 

(6 x 10 ' 4 ) ± 0.7 
20 ±1 

2.2 

(1 x 10 ~ 3 ) ± 0.7 
18 ±1 

2.5 

(4 x 10" 4 ) ±0.8 
17 ± 1 

3.0 

(4 x 10 -4 ) ± 0.9 
14.5 ±0.9 

4.2 

(4 x 10" 4 ) ± 1.2 
11.1 ± 0.7 

7.6 

c [s-*k- 3 / 2 ] 

7e,T->0 [s] 

1.8 

(4 x 10 ~ 4 ) ± 1.3 
10.4 ± 0.6 

2.0 

(2 x 10 ” 3 ) ± 1.4 
9.9 ± 0.6 

2.2 

(2 x 10 ~ 3 ) ± 1.4 
9.4 ± 0.6 

2.5 

(7 x 10 -4 ) ± 1.5 
8.7 ± 0.5 

3.2 

(3 x 10 -3 ) ± 1.8 
7.4 ±0.5 

cq 

U 

xt 

c [s” 1 k- 3 / 2 ] 

r e,T-K) [s] 

1.8 

(4 x 10 ~ 3 ) ± 2.3 

5.7 ± 0.3 

1.9 

(4 x 10 ~ 3 ) ± 2.4 
5.5 ± 0.3 

2.0 

(3 x 10 ~ 3 ) ± 2.5 
5.4 ± 0.3 

2.2 

(2 x 10 ~ 3 ) ± 2.6 
5.1 ± 0.3 

2.6 

(4 x 10 -3 ) ± 2.9 
4.6 ±0.3 

16 

xt 

c [s- 1 k- 3 / 2 ] 

Te.T-s-0 [s] 

1.8 

(5 x 10 ' 3 ) ± 2.7 
4.9 ± 0.3 

1.9 

(5 x 10~ 3 ) ± 2.8 
4.8 ±0.3 

2.0 

(8 x 10" 4 ) ±2.8 
4.7 ± 0.3 

2.1 

(4 x 10 -3 ) ± 2.9 
4.5 ±0.3 

2.5 

(6 x 10 -3 ) ± 3.2 
4.1 ±0.3 

76 

xt 

c [s^ 1 k- 3 / 2 ] 

r e,T->0 [ S ] 

1.8 

(2 x 10" 2 ) ± 13 
1.04 ±0.06 

1.9 

(2 x 10 -2 ) ± 13 
1.03 ±0.06 

1.9 

(2 x 10" 2 ) ± 13 
1.03 ±0.06 

1.9 

(2 x 10" 2 ) ± 13 
1.02 ±0.06 

2.0 

(2 x 10“ 2 ) ± 13 
1.00 ±0.06 


50 s. This would increase the density (without correc¬ 
tions) in the converter from the currently observed value, 
e ex e v (1.08 K) P(1.08 K) U CN r to t/U conv = (23 ± 3) cm" 3 . 

For instance, if r to t = 100 s were attained, a modest 
goal for the apparatus (which would make e v ~ 0.9), and 
if the CN beam were better optimized for the converter 
vessel (VcN/bconv ~ 1), then a uncorrected UCN density 
of e ex e v P rtot ~ 600 cm -3 would be expected. 

In the second part of this paper, continuous measure¬ 
ments with the converter temperature scanned from 1.2 
to 2.4 K were used to study UCN production and up- 
scattering. The three processes expected to contribute to 
the up-scattering of UCNs are: 1-phonon absorp¬ 

tion, 2-phonon scattering, and roton-phonon scattering. 

The continuous measurements for T < 1.95 K were 
fitted using the theoretical temperature-dependences of 
these processes, with both calculated values of the A, B, 
and C coefficients associated with the above processes 
respectively, as well as by varying them. Our analysis 
revealed the need to include only 2-phonon scattering, 
i.e. rjp 1 (T) = BT 7 , with a value of B = (4 16) x 
10~ 3 s _1 K~ 7 being valid in describing the data. This 
conclusion is the same as from earlier experiments [421 [601 


ESI but performed using an enlarged temperature range 
and a more thorough data-fitting procedure. Extracted 
from the analyses are the limits: C < 5.3 s _1 K -3 / 2 (95% 
C.I.), i.e. < 29 % of its theoretical value; and A < 61 s _1 
(95% C.I.), i.e. < 47 % of its theoretical value. 

The temperature-dependence of the UCN production 
rate, P re i(T), was calculated from inelastic neutron scat¬ 
tering in this paper. The calculated form of P re i(T) 
qualitatively explains the interesting kink observed in 
N C (T) at T\, and is required for obtaining a good fit of 
N C (T < 1.95 K). A discrepancy caused by an increased 
UCN loss and/or reduction in UCN production was ob¬ 
served when T > 1.95 K. This is most likely due to us not 
taking into full account the details of superfluid helium 
close to T\. 

This study demonstrates that there is still a lack of full 
understanding of the interactions between UCNs and su¬ 
perfluid 4 He. It has been suggested [53] that the observed 
lack of UCN loss due to roton-phonon scattering is a con¬ 
firmation of the proposed strong preference for co-linear 
momentum vectors in the three quasi-particle coupling 
constant when the momentum and energy transfers are 
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large (> 3.4 A 1 and > 1.1 meV respectively) [7Bj. Fur¬ 
ther theoretical work, such as that done in Ref. [77], but 
for the neutron, could shed light on the current discrep¬ 
ancy. Also, expanding on experimental work, such as im¬ 
proving the precision and temperature range of measure¬ 
ments of the up-scattered UCN spectrum [75] or under¬ 
standing the discrepancy of multi-phonon UCN produc¬ 
tion in pressurized helium with calculations 131. could 
provide valuable insights into this problem. 
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VIII. APPENDIX: TEMPERATURE 
VARIATION OF UCN PRODUCTION 

Calculations of the relative UCN production rate 
P rel (T) in Fig. [9] are performed, following Refs. m and 
in], using the expression 

P(T) = NaU c ^ f s(A n ) A n dA„ , (24) 

Jo d.A n 

where N is the helium number density, a the bound neu¬ 
tron scattering length of 4 He, U c and k c are the critical 
energy and momentum of UCN storable in the converter 
volume, A n is the incoming CN wavelength, and d(/>/dA„ 
is the incident CN particle flux spectrum. The term s(A n ) 
is defined in terms of the dynamic structure factor of 4 He 
S(q,uj), where hq and hoj are the momentum and the en¬ 
ergy transfers, via 


s(A„) = h J S(q,uj) S(hu> — E n ) du! , (25) 

where E n = /i 2 /(2m„A„) = (Hk n ) 2 /(2m„), with Hk n and 
m n being the neutron momentum and mass. 


The raw data [50] of S(q,uj) obtained from inelastic 
neutron scattering described in Ref. [73] is used in these 
calculations. These measurements were performed by 
scattering incident CN with a fixed energy of 3.8 meV 
off high-purity liquid 4 He in a 18 cm 3 container with the 
energy transfer determined using a time-of-flight neutron 
spectrometer. The measurements cover temperatures be¬ 
tween 1.24 K and 4.95 K and k between 0.3 and 2.1 A 1 . 
Plots of s(A n ) extracted from this data at selected tem¬ 


peratures are shown in Fig. 10 


The instrument energy resolution of these measure¬ 
ments is ~ 100 fj,eV (FWHM), which at the single¬ 
phonon peak, Hui ~ 1 meV, corresponds to a g-resolution 
of ~ 0.03 A . This dominates the width of the single¬ 
phonon peak because, for example, the width of the 
single-phonon peak has been measured using neutron 
triple-axis and spin-echo spectrometers to be ~ 7 /reV 
[511 . When using a broad cold neutron spectrum, as is 
done here, the effect of the energy resolution is small 
when calculating P(T). 


Linearly extrapolating s(A n ) 
and shown in Fig. [10 


0 as A n —>• 3A is done, 
This is motivated by the mea¬ 


surements of Refs. [52] and m- More sophisticated ex¬ 
trapolation schemes are available in], but at the desired 
calculation precision, they are unnecessary. 

The empirical approximation of the capture flux spec¬ 
trum d<A/dA n = (A n /1.8A) d(£/dA n from Ref. [6T] of the 
PF1B CN beam (also shown in Fig. [l0| ) is used in these 
calculations. To compute the integral of Eq. 24 numeri¬ 
cally, a constant spacing of dA n = 0.01 A was used. The 
S(q,ui) data points were 2D linear interpolated to calcu¬ 
late s(A n ). 

Using the above, as well as the helium density [83], the 
differential UCN production rate P(T, A n ) is calculated 
and also shown in Fig. 10 The familiar structure: a clear 
separation of single-phonon and multi-phonon UCN pro¬ 
duction, can be seen at the lower temperatures. As the 
temperature increases, the single-phonon peak becomes 
broader and less prominent, and disappears beyond T\. 

The P re i(T) points in Fig. [9] come from integrating 
P(T, A n ) over A n at different temperatures below 3.5 K. 
The statistical errors of the S(q,u) data [751150] are com¬ 
bined to calculate the statistical error in P re i(T). Since 


the spacing used to evaluate Eq. 24 and[25]is smaller than 
the spacing of the S(q,w) data, if a point is used multi¬ 
ple times, the errors are not combined independently (i.e. 
not added in quadrature). The statistical error in s(A n ) 
is ~ 2-10% in the regions that contribute significantly 
to UCN production. And the total statistical error in 
P re i(T) is between 0.3- 0.4%. 
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FIG. 10. (Color online) The PF1B CN capture flux d(/> c /dA n , 
s(A n ) (Eq. 251, and the differential UCN production rate 
P(T, An) used for the calculations of P le \(T), at selected tem¬ 
peratures. Each separate quantity plotted is in the same ar¬ 
bitrary units for all temperatures, but the plots for incremen¬ 
tally increasing temperatures are offset by —0.3 (arbitrary 
units) from the previous. Details of the calculations are de¬ 
scribed in the text. 
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